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EXPERIMENTAL ASSESSMENT OF THEORY FOR REFRACTION 
OF SOUND BY A SHEAR LAYER 


By 


Robert H, Schlinker and Roy K. Amiet 


SUMMARY 


A theoretical and experimental investigation was conducted to determine the 
refraction angle and amplitude changes associated with sound transmission through a 
circular, open-jet shear layer. Experiments were performed using a 0.91 m diameter 
open jet in the United Technologies Research Center (UTRC) Acoustic Research Tunnel. 
Free stream Mach number was varied from 0.1 to 0.4. 

For measurements directed toward assessment of the refraction angle correction, 
a discrete tone acoustic source with a frequency range of 1 kHz to 10 kHz, was 
situated in the airstream. Sound wave front angle change was measured for several 
axial source locations and two off-axis source positions. Refraction was observed 
to be significant at Mach numbers of 0.1 and greater. The experimental results 
were compared with an existing refraction theory which was extended in the present 
study to account for off-axis source positions. 

Good agreement between the refraction angle correction theory and experiment 
was obtained over the test Mach number, frequency and angle measurement range for 
all on-axis acoustic source locations* For off-axis source positions, good agree- 
ment was obtained at a source-to-shear layer separation distance greater than the 
jet radius. Measureable differences between theory and experiment occurred at a 
source-to-shear layer separation distance less than one jet radius. This disagree- 
ment, at present, is not viewed as a failure of the theoretical angle correction. 
Rather, it is believed due to the sensitivity of the experimental method to small 
changes in the assumed zero-thickness shear layer position at small source-to-shear 
layer separation distance. As previously predicted, experimental results indicate 
that axial variations in shear layer thickness and shear layer divergence need not 
be included in the refraction correction procedure. The frequency independence 
of the analysis was also verified. 
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Theory for the refraction amplitude change associated with sound transmission 
through a shear layer could not be verified due to limitations of the experimental 
technique. Indirect verification of the theory, however, was achieved over a 
limited range of frequency and Mach numbers using previously published experimental 
results. 

A shear layer turbulence scattering experiment was conducted at 90^ to the open 
jet axis for the same freestream Mach numbers and axial source locations used 
in the refraction study. Significant discrete tone spectrum broadening and tone 
amplitude changes were observed at open jet Mach numbers above 0.2 and at acoustic 
source frequencies greater than 5 kHz. More severe turbulence scattering was 
observed for downstream source locations. Consequently, sound radiated at angles 
close to the jet axis will encounter greater turbulence scattering than that 
radiated at 90^ to the open jet. 
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INTRODUCTION 


Open jet acoustic test facilities are currently used to investigate the 
effect of forward flight on aeroacoustic noise mechanisms. The technique permits 
studying a wide variety of problems such as airframe noise, model propeller and 
rotor noise, isolated airfoil and blown flap noise, and jet noise forward flight 
effects. For tests conducted at free stream Mach numbers less than 0.1, measure- 
ments outside the airstream can be used to directly infer the source noise 
characteristics. However, at higher Mach numbers, the open jet technique^'^ 
is susceptible to errors due to the presence of the shear layer through which the 
sound is transmitted. The shear layer serves to refract, reflect, and scatter the 
sound radiated from the model. These effects significantly alter the acoustic 
source directivity pattern and hence alter the conclusions drawn from a particular 
experiment. 

The lack of a firm understanding of these effects has created a clear need 
for a validated shear layer correction procedure. The present study was undertaken 
to experimentally assess the refraction angle change due to propagation through a 
shear layer with a lesser effort devoted to the refraction amplitude change. Exist- 
ing refraction theory was extended to treat the case of an off-axis source. Both 
on-axis and off-axis acoustic source locations were tested for a circular jet 
geometry. Acoustic source frequency and jet Mach number were also varied to 
evaluate the dependence on these parameters If the thoery could be validated, a 
correction procedure would be available for reducing far-field noise data acquired 
in open jet test facilities. 


Previous Theoretical Investigations 

Because the use of open jet wind tunnels for determining the acoustic radiation 
properties of test models is a relatively new technique, methods have only recently 
become available for correcting the acoustic data obtained in such a facility for 
the effects of refraction. However, the refraction problem had received previous 
attention because of its importance in the jet noise problem. The case of refrac- 
tion of a plane sound wave by a plane, zero thickness shear layer was first 
correctly treated by Ribner (ref. 1) and Miles (ref. 2). Gottlieb (ref, 3) extended 
the analysis to the case of a point source beneath a plane shear layer. Graham 
and Graham (ref. 4) considered the problem of a plane wave interacting with a finite 
thickness shear layer and subsequent publications of theirs consider the field of 
specific singularities near a shear layer. Amiet (ref. 5) and others also considered 
the sound field of specific source types near a shear layer. 
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The above studies were not concerned with correcting acoustic wind tunnel data, 
however, and further analysis of the problem was needed. In particular, a method 
was required which made no assumptions regarding the nature of the sources (e,g., 
monopole, dipole, etc.), since, in general, the nature and distribution ofl 
sources in a given test are unknown. Subsequently, Amiet (ref. 6) derived a correc- 
tion procedure which applied both an angle change and an amplitude change to the 
data. This assumed a plane zero thickness shear layer. It was only by this 
separation of angle and amplitude effects that it was possible to arrive at a 
correction independent of source type. The problem was also analyzed by Jacques 
(ref. 7) who arrived at the same result for refraction by a plane shear layer, and 
also considered the case of a source on the centerline of a cylindrical shear layer 
of zero thickness. The solution for refraction by a thick cylindrical shear layer 
with a source on the centerline was given by Tester and Morfey (ref. 8) and for a 
thick plane shear layer by Amiet (ref. 9). This work of Amiet also gives a review 
of the several correction procedures. 

The above corrections are in algebraic closed form. Candel (ref, 10) recently 
developed a numerical ray tracing procedure which appears to give results very close 
to the closed form solutions. Recent studies by Tester and Burrin (ref, 11) 
indicate that the axial variation of the shear layer has little effect on the 
refraction corrections. Thus, it is possible to model the open jet shear layer 
without axial variations. 

References 6-11, above, rely on the technique of correcting both angle and 
amplitude, Mani (ref. 12) employed a different approach to the problem, but this 
appears to be less well established than the above mentioned techniques of correcting 
angle and amplitude. 


Previous Experimental Investigations 

Few experimental studies exist documenting the shear layer refraction and 
scattering characteristics. Amiet reported the results of an exploratory 

experiment for the angle correction which tended to support the theoretical pre- 
dictions in reference 6. Candel et al. (ref. 13) provided an experimental assess- 
ment of a refraction theory in addition to presenting preliminary results for 
turbulence scattering in an open jet. Similarly, Ahuja et al. (ref. 14) recently 
reported an assessment of the refraction theory of reference 11. Also included was 
an exploratory study of turbulence scattering at low jet Mach numbers and low 
acoustic source frequencies. All of the above experimental investigations were 
limited to on-axis source locations. 
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Method of Approach 


In the present study, the refraction angle change was assessed by cross- 
correlating far-field microphones to measure the local sound wavefront propagation 
angle outside the open jet. Knowing the propagation angle, the acoustic ray 
normal to the wavefront was traced back to the shear layer crossing point. In 
determining this point, the shear layer was assumed to be cylindrical with a radius 
equal to the inlet nozzle radius and to have zero thickness. From this, the 
original propagation angle inside the shear layer was determined for comparison 
with refraction angle correction theory. The original analysis of Amiet presented 
in references 6 and 9 was extended here to accommodate off-axis acoustic source 
locations. The refraction amplitude correction was assessed with a dipole acoustic 
source located in the open jet airstream. The difference between the measured 
far-field tone directivity and the theoretical dipole directivity pattern inside 
the jet represents the amplitude correction. Here, the theoretical angle correction 
was used to link propagation angles inside and outside the airstream. The experi- 
mental amplitude corrections were then compared to the refraction theory discussed 
later. 

A simple experiment was also conducted to assess the importance of turbulence 
scattering. A discrete tone was used to investigate spectrum broadening effects 
as well as changes in the tone amplitude. Various frequencies and Mach numbers were 
employed to determine when scattering effects become significant. The importance 
of turbulence length scale changes on the scattering phenomenon was also investi- 
gated by locating the acoustic source at various axial stations. 
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Acoustic source frequency, Hz 
Microphone sensitivity 
1- h/a 

Cross-correlation amplitude 

Source-to-shear layer distance. Also corresponds to source- to- 
lip-line distance. 

Distance between source and quarter velocity point in shear layer 
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Unit vector from retarded source position to shear layer crossing 
point 
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DESCRIPTION OF THE EXPERIMENT 


Acoustic Research Tunnel 

The experimental study was conducted in the UTRC Acoustic Research Tunnel. 

A detailed description of the facility is given in reference 15. The tunnel, 
shown schematically in figure 1 is an open-circuit, open- jet design. The inlet 
is provided with a high length-to-diameter ratio honeycomb section and a series 
of turbulence suppression screens. These features, in conjunction with a 
large tunnel contraction, provide a spatially uniform, temporally steady flow 
with a controlled test section turbulence level of approximately 0.2 percent. 

An inlet nozzle with a radius, R^ = 0.455 m, was employed in the present test 
program. Turbulence generators and grids can be inserted upstream of the nozzle 
to provide various wake profiles (ref. 16) and a range of turbulence levels 
(ref. 17) in the test section. 

The open jet test section is surrounded by a sealed anechoic chamber 4.9 m 
high, 5.5 m long (axial direction,), and 6.7 m wide. The interior walls of the 
chamber, shown in figure 2, are lined with 0.3 m fiberglass wedges. The chamber 
has been found to be anechoic (ref. 15) for broadband noise over a 200 Hz to 
20 kHz range of calibration frequencies. That is, the sound pressure followed a 
6 dB decay curve per doubling of source-observer separation distance within approxi- 
mately 1/2 dB over this frequency range. 

Since pure tone source frequencies were used in the present experiment, 
intensity versus distance calibration measurements were made to verify the free 
field behavior. For these tests a speaker was located on the centerline of the 
open jet test section. A microphone was traversed radially outwards from the 
acoustic source at various angles to the jet axis. Figure 3 shows a typical 
result at 90^^ to the jet centerline. The intensity changes followed the 
inverse-square law decreasing approximately 6 dB with doubling of distance over 
the 1 kHz to 10 kHz frequency range. Similar results were obtained at 
f = 12.5, 16, 20 and 25 kHz. These tests verified the anechoic characteristics 
of the facility for pure tones over the frequency range employed in the present 
study. It should be noted that a pure tone provides a more severe test for 
evaluating free field characteristics than random noise. 

The test section airflow enters the diffuser by a collector that has 
anechoic treatment on its flow impingement lip. Initial facility tests 
reported in reference 15 identified an acoustic coupling between the inlet 
nozzle and the collector lip resulting in edge tones at high tunnel speeds. To 
suppress this noise mechanism, triangular tabs (see figure 4) were distributed 
around the nozzle periphery to disturb the azimuthal symmetry of the shear 
layer and prevent the generation of feedback tones. On the other hand, the 
shear layer thickness is increased by the tabs as will be discussed in the 
section on Shear Layer Growth. Despite this alteration of shear layer width, 
the measured angle change associated with the wavefront propagation through 
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the finite thickness shear layer showed good agreement with the change calculated 
using the circular vortex sheet theory developed in the Resent study. 

The diffuser in figure 1 operates unstalled and is thus not a major source 
of background noise. To avoid tunnel fan noise from propagating upstream into the 
anechoic chamber a Z-shaped muffling section with two right angle bends and 
parallel treated baffles is located between the diffuser and the fan. The 1500 hp 
centrifugal fan exhausts to the atmosphere through an exhaust tower. 

Tunnel speed is determined from total pressure measurements at the contraction 
inlet and static pressure measurements within the sealed anechoic chamber. Since 
losses are confined to the boundary layer, total pressure upstream and downstream 
of the contraction are predicted and have been verified, to be equal. The test 
section velocity has been shown to be temporally steady. 

Experimental Arrangement 

Refraction Angle Correction Experiment - Figures 2 and 4 show the anechoic 
chamber test arrangement employed for the angle correction measurements. The 
nozzle inlet and collector described earlier are evident in the photos. 

An acoustic source, consisting of a speaker enclosed by an aerodynamic 
fairing (see f^ig. 4), was located in the open jet airstream. The fairing measure- 
ments were: 17 cm long (axial direction), 10 cm high, and 5 cm deep. Leading and 

trailing edges of the fairing were rounded. The wake shed by the body produced 
no significant discrete vortex shedding tones in the far-field acoustic spectrum. 
Broadband wake noise did not influence the present experiment which employed only 
pure tones. 

The speaker fairing was held in place by a support arm which extended through 
the open jet to a rigid stand outside the flow. The support arm consisted of a 
short section of cylindrical tubing and an aerodynamically shaped fairing. 

Although the cylinder generated an aeolian tone, the dipole radiation pattern was in 
a vertical direction perpendicular to the horizontal plane containing the speaker 
and the far-field microphones. Therefore this extraneous noise mechanism was not 
sensed. Similarly, noise generated by impingement of the highly turbulent shear 
layer on the aerodynamic fairing radiated in a vertical direction and was not 
sensed by the far-field microphones. Acoustic insulation was applied to the support 
stand and the exterior of the inlet nozzle. 

The speaker face, with its 2.5 cm diameter aperture, was covered with a plastic 
grid of the same aperture size. The grid prevented flow separation and distortion 
of the mean velocity field at the speaker-flow interface. The center of the grid 
aperture marked the acoustic source position. All far-field microphone locations 
were referenced to this station. 
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The acoustic wavefront angle outside the open jet was measured by! a pair 
of far-field microphones. The microphones, designated as m , and m in figures 

i 2 

4 and 5, were mounted at grazing incidence on a rotating boom. The boom pivot 
point, located beneath the open jet, was aligned with the acoustic source 
position. Although figure 5 shows the source situated on the tunnel centerline, 
of f~centerline stations were also investigated. The microphones scanned a 
horizontal plane at the same height as the tunnel centerline and the acoustic 
source. Different microphone separation distances, iy were used depending on 
the acoustic source frequency. The fixed microphone radius, r = 1.83 m, was 
measured to the midpoint of the microphone array. The radius was selected to 
permit scanning the boom close to the jet axis while still maintaining a large 
separation distance between the microphones and the downstream collector. This 
avoided potential near field reverberation problems when the acoustic source 
was situated at the furthest downstream location. The microphones were mounted 
rigidly using guy wires and braces to minimize any vibration or deflection when 
the open jet was operated. 


The boom angular position was defined by the angle, 0^, measured from the 
downstream jet axis. The angle was referenced to the midpoint between the two 
microphones. A line connecting the microphones formed approximately a 90° 
angle with the radius, r, A potentiometer, attached to the boom pivot shaft, 
provided a voltage output corresponding to the boom angular position. The 
voltage was amplified and displayed on a digital voltage readout providing a 
resolution of 0.016 deg. The system was calibrated daily using the open jet 
centerline and a line perpendicular to the centerline as references. The 
accuracy of the calibration was approximately 0.3 degrees. 

Refraction Amplitude Correction Experiment - To provide a dipole noise 
source, a segment of polished cylinder was mounted transverse to the open jet. 
The 1.27 cm diameter cylinder, shown schematically in figure 6, was centered 
vertically on the jet axis. The cylinder was held in place by a larger 2.54 cm 
diameter, cylinder which extended through the open jet shear layer to an 
external support stand. The lower end of the test cylinder was held rigid by 
several guy wires connected to points beneath the open jet. None of the wires 
crossed the flow field upstream of the cylinder. 


A splitter plate was attached to the rear stagnation point of the support 
cylinder to suppress its vortex shedding mechanism. The guy wire vortex 
shedding frequency was approximately a factor of 10 higher than the test cylinder 
frequency. Consequently, the dominant tone in the acoustic spectrum was provided 
by the test cylinder. 


The cylinder directivity pattern outside the flow was measured using the 
previously described boom arrangement shown in figure 5. The pivot point was 
aligned with the test cylinder centerline. Microphone m 2 was used to measure 
the sound pressure level as a function of measurement angle, 6^, This angle 
was now referenced directly to m^ instead of the midpoint between the microphones. 


J '' 
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Sound generated at the, cylinder ends propagated through the shear layer at 
almost the same angle as sound generated at the center of the cylinder. This 
is illustrated in figure 6 with microphone m situated at 90^ to the jet axis. 

The acoustic ray reaching m 2 from the lower end of the cylinder propagates 
through the shear layer at 12.5 degrees relative to the local normal. 

Discrete Tone Turbulence Scattering Experiment - An acoustic source, 
consisting of a speaker with a 9 cm diameter exponential horn, was located on 
one side of the open jet (fig. 7). Discrete tone sound waves were transmitted 
through both open jet turbulent shear layers to an array of five microphones orien- 
tated at grazing incidence on the opposite side of the airstream. The transmitted 
tone sound pressure level and spectrum shape were monitored as a function of 
increasing tunnel Mach number at various axial stations, X. The speaker and 
microphone array formed a horizontal plane which contained the open jet i 
axis. The speaker centerline formed a 90^ angle with the jet axis. The 
microphone separation distance, £, measured relative to m in figure 7, was 
determined by the jet Mach number. This will be discussed later. 


Instrumentation 

The various speakers were powered by a 150 watt single channel amplifier. 

A sine wave generator provided the amplifier input signal. The generator 
frequency was held to within 2 Hz of each source frequency evaluated in the 
test program. 

Far field acoustic sound pressure levels were measured with commerically 
available 0.635 cm diameter condenser microphones at grazingj incidence. Since 
microphone m 2 penetrated the shear layer in certain test configurations, 
protective grids were used at all times. Admittedly, the microphone frequency 
response was not flat due to the protective grid, but, comparisons of absolute 
sound pressure level measurements at different frequencies were not needed in 
the present investigation. For example, the angle correction experiment 
measured the acoustic wavefront propagation angle which is independent of 
the radiated sound pressure level. The amplitude correction experiment used 
only the normalized acoustic source directivity pattern at each frequency. 
Finally, the scattering experiment monitored the change in acoustic sound 
pressure level for a fixed discrete tone as the open jet Mach number increased. 
Absolute sound pressure level comparisons between different frequencies were, 
therefore, not needed. Based on the above arguments atmospheric attenuation 
was not relevant to the present study. All microphones were calibrated daily 
with a 250 Hz pistonphone. 

Microphone signals and the sine wave generator signal were recorded on 
magnetic tape during certain portions of the test. The frequency response of 
the FM tape system was flat over the frequency range investigated here. 
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Spectrum analysis was conducted with a narrow bandwidth, 500 line, real 
time spectrum analyzer-ensemble averager. Each spectrum was generated from 
256 ensemble averages. The effective noise bandwidth of the analyzer filter 
influenced the discrete tone scattering results, the consequences of which will 
be discussed later. 

Cross-correlations were performed with a 400 line, real time correlation 
and probability analyzer. Input signals were not filtered. 

Mean velocity measurements of the open jet shear layer were acquired with 
a 0.025 mm diameter hot film operating in conjunction with a constant 
temperature anemometer system. The probe was calibrated in the tunnel test 
section. 


Test Program 

Refraction Angle Correction Measurements - Measurements were performed at 
various acoustic source locations, open jet Mach numbers, and source frequencies 
to assess these parameters in the refraction angle correction theory. Two 
coordinates defined the source location: the axial station, X, and the source- 

to-open jet-lip line separation distance, h. Both coordinates are shown in 
figure 5. Table I: lists the test conditions evaluated in the present study 
with the source position coordinate normalized by the open jet radius. In- 
cluded in the table is the ratio, h/r, which defines the far-field microphone 
position relative to the source-to-lip line distance. This parameter is an input 
to the experimental and theoretical angle and amplitude correction calculations. 

The lip line also defines the location of the zero thickness shear layer used in 
the refraction theory described in the section titled "Theoretical Formulation 
of the Refraction Problem". Therefore, h is at times referred to as the source-to- 
shear layer separation distance. It should be noted that since the source-to- 
microphone distance, r, was fixed, changes in h resulted in both h/R and h/r 
changing . 

With the acoustic source situated at X/R = 1.33 and h/R = 1, the measure- 
ment angle, 0^, was varied from 22.5 to 105° in 7.5^ increments. Angles 
less than 22,5° were not investigated since microphone m^ in figure 5 would 
then penetrate the open jet shear layer. On the other hand, angles greater than 
105° were not investigated since an adequate separation distance between m^ and 
the upstream anechoic chamber wall must be maintained. Line-of-sight inter- 
ference by the tunnel inlet between the source and microphone m does not 
occur until 0^ exceeds 136°. Moving the source to X/R = 2.66 and 4.0 permitted 
increasing the upstream measurement angle to 127.5°. However, downstream 
measurement angles less than 30° could not be attained due to shear layer 
spreading. 
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The same measurement angle range was used for the of f-centerline source 
positions listed in Table iL For h/r =0.36, measurement angles less than 30^ 
could not be investigated since m 2 penetrated the flow. As in the on-centerline 
test geometries, line-o£-sight interference was not a problem. 

Refraction Amplitude Correction Measurements - The far-field directivity 
pattern of a dipole source was measured at two source locations and two open 
jet Mach numbers to assess these parameters in the refraction amplitude 
correction theory. Changing the Mach number also permitted evaluating the j- 
acoustic source frequency dependence. The test geometry, listed in Table II, 
coincided with two of the test configurations used in the angle correction 
experiments. Identical conditions were needed since the angle correction must 
be applied to the dipole directivity before the amplitude correction can be 
isolated. The dipole directivity was measured at a minimum of twelve angles 
over the range 0^ = 25° to 110°. 

Note that the off-center source location corresponding to Test Configuration 
7 in Table Il^represents a source displacement towards the observer. This dis- 
placement direction was selected based on the calculated amplitude corrections 
obtained from the theoretical analysis to be discussed later. It was determined 
that the amplitude correction is the largest as the source approaches the shear 
layer closest to the observer. This was expected to facilitate the experimental 
assessment of the amplitude correction. 

Discrete Tone Turbulence Scattering Measurements - Discrete tone scattering 
effects due to shear layer turbulence were evaluated at the Mach numbers 
investigated in the angle correction experiment. The tone frequencies corres- 
ponded to f = 1, 2.5, 5, 10, 17.5, and 25 kHz, The importance of changes in 
turbulence length scale was also investigated by locating the source at three 
axial stations given by X/R^ = 1.33, 2,66 and 4,0 
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DEFINITION OF THE OPEN JET SHEAR LAYER 


The objective of the mean velocity measurements described here was to document 
the shear layer characteristics of the UTRC Acoustic Tunnel as a function of axial 
position and Mach number. These measurements were obtained using hot wire anemometry. 
Radial distributions of the mean velocity will be presented first to illustrate how 
the shear layer features change with axial distance. Using these results, the 
shear layer growth rate will be presented as a function of axial location. Finally, 
to isolate changes in the velocity distribution with axial location, the velocity 
profiles are plotted in a dimensionless form using a similarity parameter based 
on the shear layer growth rate. The resulting curves from various axial stations 
collapse onto one curve indicating the profiles are similar. This similarity, which 
has been well substantiated by other investigators, permits comparisons of the UTRC 
open jet test facility with other facilities. 


Mean Velocities 

Mean velocity measurements were obtained at an azimuthal position corresponding 

to the position of one of the tabs on the tunnel inlet nozzle. Figure 8 shows the 

radial distribution of the mean axial velocity at various stations downstream of the 

nozzle exit. Here, the local velocity, U, is normalized by the open jet velocity, 

U , while the axial position, X, and radial position, R, are normalized by the jet 

radius, R . The four Mach numbers, M = 0.1, .2, .3 and .4, in each plot, correspond 
o 

to the test conditions used in the refraction and scattering studies. 

The potential core and mixing zone are clearly identified. The distributions 
exhibit a familiar behavior with the velocity profiles spreading radially with 
increasing axial distance. The bulge in the velocity profiles at X/R = 0.66 and 
1.33 was due to the severe mixing forced by the tabs attached to the open jet 
nozzle. This feature is smoothed out with increasing distance from the nozzle exit 
disappearing altogether beyond X/R =1.33. A comparison of the curves at a fixed 
value of X/R indicates that the normalized velocity profiles are independent of 
Mach number at each axial station. 


Shear Layer Characteristics 

Figure 9 shows a schematic representation of the coordinates needed to describe 
the shear layer characteristics. Since there is no universally accepted definition 
of shear layer boundaries, the inner and outer edges of the shear layer are 
defined here by the 90 percent (U/U = 0.9) and 10 percent (U/U = 0.1) velocity 
lines, respectively. The distance ?rom the jet centerline to tSe half -velocity 
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line (U/U^ = 0.5) is denoted by %/2* local shear layer thickness, 5, is 

defined as the difference between the inner and outer edges of the open jet. 

The average radial position of the inner and outer edge of the shear layer 
was determined from the velocity profiles in figure 8, The result is shown in 
figure 10 in addition to the average half -velocity position. The shear layer 
thickness was then calculated as a function of axial position using these average 
characteristics. This result is shown as open circles in figure 11 in addition 
to an empirical straight line fitted through the data points downstream of X/Rq = 
1.33. The straight line could not be extended upstream since it would predict 
a shear layer thickness which exceeds the value measured at X/R^ = 0.66. 

Also shown in figure 10 is the well-known analytical solution by Goertler (see 
page 96, ref. 18) describing the growth of the shear layer in a two-dimensional 
plane jet. The analytical solution was used to calculate the inner and outer 
edges of the shear layer based on the 90 percent and 10 percent velocity lines. 

The Goertler solution is adequate for determining the growth rate of axisymmetric 
jets. This is evident in the good agreement between theory and the measured open 
jet shear layer thickness (see figure 11) reported by Gandel, et al. (ref. 19). 

Note that tabs were not necessary in the facility described in reference 19. 

Several conclusions can now be drawn regarding the UTRC open jet. First, 
the shear layer thickness close to the nozzle exit is large when compared to other 
facilities. This is due to the early mixing forced by the tabsi attached to the 
nozzle exit. Second, the shear layer thickness grows at a slower rate when compared 
to the standard jet without tabs. At large X/R , shear layer thickness approaches 
the Goertler solution. Despite the alteration of the UTRC shear layer by the tabs, 
the acoustic refraction and scattering measurements showed good agreement with 
studies conducted in other facilities. This will be discussed later. 


Similarity of Velocity Profiles 

To permit direct comparison of the open jet velocity profiles at various 
axial stations, the mean velocity distributions in figure 8 were plotted versus 
the similarity parameter, n. Here t> is defined as, p = (R - R]^/2)/<5. A comparison 
of the normalized velocity profiles is given in figure 12 for three different 
axial stations at Mach number, M = 0.4. The various curves at the three stations 
collapse onto a single curve. Also shown is the standard hyperbolic tangent 
analytic curve used for two-dimensional plane-jet shear layers. The good agreement 
between experimental data and the hyperbolic tangent curve means this relationship 
adequately' descri.bes the UTRC shear layer characteristics. This relation can 
therefore be used by other investigators in comparing their shear layer properties 
with those of the UTRC Acoustic Tunnel. 
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THEORETICAL FORMULATION OF THE 
REFRACTION PROBLEM 


Review 


The basic approach for most of the present shear layer correction theories 
was initially given by Amiet (ref. 6). This analysis treated the case of a 
source beneath a plane zero-thickness shear layer. The observer was assumed to 
be in the geometric and acoustic far field of the source. However, the shear 
layer could be at any distance from the source, which needed be neither acoustically 
compact nor a point source. This analysis, based on the previous work of Ribner 
(ref. 1) and Miles (ref. 2), requires that the measurements made outside the 
shear layer be corrected in both angle and amplitude. The result is a correction 
method which is independent of frequency and source type (i.e., monopole, dipole, 
quadrupole) . 

The same techniques were applied by Jacques (ref. 7) to the case of a 
source on the centerline of a cylindrical shear layer. This case, although 
quite similar in principle, does require additional assumptions. In particular, 
it must be assumed that there are no multiple reflections within the shear 
layer; that is, the waves reflected from the shear layer are assumed to be 
completely damped before again reaching the shear layer. Also, the shear layer 
must be in the acoustic far field of the source; if the shear layer were not in 
the acoustic far field, then as shown by Amiet (ref. 9), the shear layer correction 
would depend on source type and frequency. That is, the solution of Amiet is not 
valid at low frequency where the wavelength of sound is comparable to the jet 
diameter. For open jet diameters encountered in practice, this is not usually a 
severe limitation. The cutoff frequency (frequency below which the surrounding 
chamber becomes non-anechoic) of open jet acoustic test facilities typically 
correspond to a wavelength to diameter ratio on the order of 1. 

The high frequency case of a source on the centerline of a cylindrical shear 
layer with a thickness significantly greater than a wavelength was treated by 
Tester and Morfey (ref. 8), Although not explicitly derived as a correction for 
open-jet wind tunnel measurements, the results of the analysis can be applied 
toward this end. 

It is interesting to note that the shear layer corrections for the two 
cases of zero thickness and infinite thickness shear layers differ by only a 
few hundredths of a decible over most of the angular range as shown in reference 
9. This would be expected since the basic difference between the two corrections 
is that the zero thickness case assumes a reflected wave exists while the infinite 
thickness case assumes that all the energy is transmitted. Since the amplitude 
of the reflected wave is small over most of the angular range (see e.g., reference 
6) and since the reflected energy is proportional to the square of the amplitude. 
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one should expect little difference between the corrections. At the extreme 
angles where the reflected wave can have a significant magnitude, both the zero 
thickness and the infinite thickness corrections should be used with caution since 
neither correction accounts for multiple reflections within the shear layer. Thus, 
the correction procedure for an infinite thickness shear layer used by Ahuja et al 
(ref. 14) has the same angle correction as that of Amiet (ref. 6 and 9) and 
Jacques (ref. 7) and an amplitude correction differing by a few hundredths of a 
dB over most of the angular range. It is, for all practical purposes, equivalent 
to the correction procedure for a zero thickness shear layer. A more detailed 
comparison between the various zero thickness, infinite thickness, plane shear 
layer and cylindrical shear layer results is given in reference 9. 

Shear Layer Correction for An Off-Axis 
Acoustic Source 

For the case of a source on the centerline of a cylindrical shear layer, as 
was discussed in reference 9, the shear layer must be in the acoustic far field 
of the source for the shear layer correction to be independent of source type 
and frequency. Thus, ray acoustics principles are applicable and the shear 
layer correction can be derived by ray tracing arguments. For the case of a 
source off the centerline, the derivation will also be based on ray acoustics. 

The same far-field assumption required for the on-axis source solution applies 
here. The main difference between this case and the on-centerline case is that 
the algebra for the ray tracing of the off-center case becomes significantly more 
complex. The basic principles governing the crossing of the shear layer by a 
ray remain the same, however. 

Figures 13 and 14 show the basic geometry and the coordinate system used in 
the derivation. The zero thickness cylindrical shear layer has a radius, a, 
and the source is located off the centerline by a distance a-h. A Cartesian 
coordinate system has the origin located at the source, S, the x axis parallel 
to the tunnel axis and the z axis along the normal from the tunnel centerline to 
the source. A ray travels from the source along the path labeled r-j^ crossing 
the shear layer at point x^ , continues along path r^ to the 

observer 0 at x-, ^2^ * ^ 2 ^ * plane tangent to the shear layer at point 

X. , y * , z ’ makes an angle a with the x-y’ plane. Angles 0^and c|)^ describe 
tfte ray paih beneath the shear layer, 6|.and cj)^ describe the ray path above the 
shear layer and 6^ and (j>^ describe the observer position relative to the source. 

The angle 0 is the angle between a given line and the x axis or a line 
parallel to the x axis. Thus 0^ is the angle between the ray beneath the shear 
layer and the x axis, 0^ is the angle between the transmitted ray and a line 
parallel to the x axis at the shear layer crossing point, and 0^ is the angle 
between the x axis and the line joining source and observer. The angle cj) is the 
angle between the x-y plane and a plane containing a given line and parallel to 
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the X axis. Thus, is the angle between the x-y plane and the plane defined by 
the X axis and the ray beneath the shear layer, (j)^ is the angle between the x-y 
plane and the plane containing the transmitted ray and parallel to the x axis, 
and c() is the angle between the x-y plane and the plane containing the x axis 
and tSe observer point. The unprimed coordinate system is defined by rotating the 
X, y’ , coordinate system through an angle a about the x axis so that z is 
normal to the tangent plane. 

The detailed derivation of the off-axis source correction is given in 
Appendix A. Results of this calculation can be summarized as follows; the 
amplitude correction is different from the on-axis case for all source-observer 
positions. In the special case of an observer located in the plane containing 
the source and open-jet centerline, the angle correction is identical to that 
obtained for an on-axis source or for the case of a plane shear layer. 
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EXPERIMENTAL ASSESSMENT OF REFRACTION THEORY 


Angle Correction 

Experimental Approach - The primary objective of the experimental program was 
to measure the angle change associated with acoustic wave-front propagation through 
an open jet shear layer. While various combinations of source-observer locations 
could be evaluated the present study was limited to acoustic source and far-field 
microphone positions situated in a plane which coincided with the jet axis. Under 
this conditions, angles = 90 and a = 0 in figures 13 and 14. The resulting 
propagation path, SBO in figure 14, was thus confined to a plane normal to the 
tangent plane at the shear layer crossing point. 

The three-dimensional geometry in figures 13 and 14 can therefore be replaced 

by the two-dimensional geometry in figure 15. Here wave fronts propagating in the 

direction, 0 , inside the air stream propagate in the direction, 9^, outside the 

open jet. T^e actual path of an acoustic ray is described by the points SBO. The 

change from 0 to 0^ is a result of refraction by the shear layer. Without flow 

the sound propagates at angle 0 . 

m 

The experimental refraction angle correction for the geometry in figure 15 was 
obtained using the test set-up shown in figure 5. Acoustic ray angle changes were 
determined from measurements of the constant phase or wave-front surfaces outside 
the airstream. In this region, the acoustic ray is perpendicular to the wave-front. 
Once the wave-front angle at the observer location is known, the ray normal to the 
wave-front ray can be traced back to the point, X , at which the sound emerged from 
the shear layer. The original propagation angle xnside the open jet is then 
obtained from the relationship 0^ = arctan (h/X^) . 

The wave-front angle was measured in the following manner. Sound was generated 
at discrete frequencies by an acoustic source inside the open jet. The phase 
difference, between the signals arriving at microphones m- and m^ in figure 5, was 
then determined by cross-correlating the signals. From this, the wave-front angle 
relative to the microphone array was calculated permitting the wavenormal ray to be 
traced back to the shear layer crossing the point. 

Measured Phase Difference and Calculated Angle Change - It now remains to derive 
the equations linking the measured phase difference and the original propagation 
angle, 0 , inside the free jet. The modeling problem is shown in figure 16. For 
M = 0, tfee angle between the wave front arriving at m 2 and the line connecting the 
array is denoted by y . This parameter accounts for small differences in the dis- 
tance, Ar, between the source and the individual microphones. This difference is 
due to inaccuracies in placing the two microphones at exactly the same radial dis- 
tance from the source. Typically, the difference was less than 0.25 cm. Assuming 
this value of Ar yields a maximum value of y^ = 0.016 radians for the smallest 
separation distance, £, shown in figure 16. 


21 


If the microphone array skewness were neglected, errors would arise in the 
phase shift measurement. For example, if Ar = 0.25 cm, then at 10 kHz a 0.07 cycle 
phase shift occurs between and m^. Since the present experiment required an 

accuracy of 0.02 cycles it was necessary to account for the skewed microphone 
alignment. That is, the measured phase shift at M = 0 was subtracted from the phase 
shift obtained at finite Mach number. This effectively calibrated the measurement 
of phase shift to account for microphone placement error. The mathematical details 
are presented below. 

The phase difference (expressed in cycles) between constant phase points 
arriving at m^ and m^ is given by the relation 


A-r= A+(A,- Aaf 


( 1 ) 


fi 

A= -3^ sin/io 


( 2 ) 


The parameters, A|^ and account for the phase lag between the acoustic input 
signal to each microphone diaphragm and the final voltage signal applied at the 
correlator. Nominally, this phase lag is less than 0.1 cycles for the frequency 
range employed in the present angle correction experiment. The microphone cartridge 
phase lag characteristics are included in this parameter. 

The parameter, ZL, was measured directly by cross-correlating the periodic 
signals arriving at the adjacent microphones. Figure 17 illustrates the correlator 
time delay trace. Based on the phase shift of the first maximum, = T*f. The 
magnitude of t has been exaggerated to make it perceptible in the figure. 

Nominally, ZL is less than 0.1 cycles. The location of the first maximum in the 
correlation function was defined as the mid-point between the zero cross-over points. 

In figure 17 the mid-point of the second cycle and a knowledge of the period were 
used to locate the mid-point of the first cycle. By using two cross-over points 
irregularities in the correlation function were averaged. The term "irregularities” 
is used here to describe a deviation of the measured correlation functions from a smooth 
sinusoidal behavior. These deviations occurred because both microphones sensed 
open jet shear layer background noise as well as source noise. Such irregularities 
in the correlation function only occured, therefore, at non-zero Mach numbers. 

For M ^ 0, the phase difference, therefore, is given by 

a \= a'+ (Aj-Ag) 

where 

= sin(/Xo + /A,) 


(3) 


(4) 
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Here describes the angle between the two wavenormal rays in figure 16. Solving 
for the difference - Aj, gives 


a't ~ At = a' - A = 



sin(/io +/!,)- sin 


(5) 


Appendix B shows how this expression may be simplified further if y is small. 
Solving for y^^ gives 

Ml = sin“'[“ff^ (a'j-At)] (6) 


From this 0^ can be determined using the relationship 


tan"*(h/Xo) 


(7) 


where 


Xq= r cos 6 ^ - (r sin h ) cot ( Mi'*’ ^ 


( 8 ) 


The expression for is derived in Appendix C, 

The above approach was successful only at M = O.l, For higher open jet Mach 
numbers, the facility random background noise dominated the cross-correlation func- 
tion as is illustrated by the two-microphone correlation function in figure 18. 
Although a periodic waveform is present in the cross-correlation function, the 
measured result is not suitable for solving equation 6 for the wavefront propagation 
angle. To verify that the time delay trace is dominated by facility noise, a second 
trace is shown with the acoustic source turned off. 

The background noise dominated the cross-correlation function despite the 
12 dB signal to noise ratio shown in figure 19. This occurred because the micro- 
phone input signals to the correlator were not filtered permitting the complete 
background noise signal to contribute to the cross-correlation. It may appear 
that this problem can be circumvented by narrowband filtering the microphone spec- 
trum shown in figure 19 to pass only the pure-tone component. However, this still 
retains in the correlation function a contribution from the cross-correlation of 
the filtered narrowband microphone signals. The periodicity of this contribution 
corresponds to the filter cutoff frequencies . Since these frequencies are close to 
the acoustic source frequency, a portion of the narrowband random contribution 
would remain in the cross-correlation. 
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An alternate cross-correlation technique developed by Schlinker (ref. 21 ) 
was therefore employed. The method involves cross-correlating the individual 
microphone signals with the periodic input signal to the acoustic driver system. 

The random signals due to facility noise are rejected leaving only the periodic 
component. The success of the method is illustrated by the correlation traces in 
figure 20. Here and 82^3 represent the individual microphone-signal 

generator cross-correlation functions. To verify that the traces are independent 
of the facility noise, one trace is shown with the source turned off. 

The mathematical details linking the correlations in figure 20 with the micro- 
phone orientation angle + y are presented in Appendix D. The final result for 
the phase difference at non-zero Mach number is 


^T=Vr’72 


(9) 


where n* represents the phase difference between the acoustic driver signal and 
the microphones. Equation (9), when used in conjunction with equation ( 6 ), gives 

^i=sin'[-^(i7;-77'2-AT)] 1 do) 

A brief discussion of the measured phase difference is now warranted since the 
angle correction is derived from this parameter. Figure 21 shows the experimentally 
determined phase difference, for a 2.5 kHz acoustic source frequency with 

the open jet operating at M = 0.3. The source position corresponded to Test 
Configuration 1 in Table I. Two different microphone spacings were used in the 
measurement. For angles near 90 , where refraction is minimal, a large separation 
(£ = 0.48 m) was used to obtain a measurable phase shift. In contract, at angles 
close to the jet axis where the phase shift is large a smaller separation (£ = 0.3 m) 
was used. 

The measured data in figure 21 are compared to the phase difference calculated 
from the refraction theory developed in the previous section. The figure shows 
excellent agreement between the theory and experiment over the complete range of 
measurement angle. The measured phase was then used to calculate the angle correc- 
tion associated with the transmission of sound through the shear layer. 

Plane Wave Propagation Assumptions - The above experimental approach assumes 
that the wave front arriving at the far-field microphone array is plane. The 
possible error in this assumption was evaluated by calculating the time for the 
sound to travel from the source to the individual far-field microphones. The 
difference between these times was then compared to the time calculated using the 
plane wave assumption. 
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For a point source in a flow, the phase, 0% of the waves is 


0 * t + 




[mx- 


>/; 


2 + 




')] 


(O 


( 11 ) 


If 0| = 0 at t = 0, then the same wave front will reach the point (X , 0, h) at 

o 

time t where 
o 




+ MX, 


=/; 


2 + 




( 12 ) 


Here is the shear layer crossover in figure 15. Using tan 9^ = h/X^ along with 
equation A. 18 relating 9^ and 9^ gives 


h(l-MCOS^t) 

Co^t 


(13) 


The time for the wave to travel from Xq to the observer at 0 (see figure 15) is 


rsin^m-h (14) 

casing; 

The total time for a sound wave to travel from the source to a far-field 
microphone is then ^2 ~ ^1’ exact time difference, ^t^, for the sound 

to reach the two separate microphones is obtained from the difference between the 
value of t^ calculated for each microphone. 

On the other hand, the approximate calculation treating the wave front as 
plane gives for the time difference (see figure 16) . 

Al„ = sin(e, -e„,, = -sj-siriAij, (15) 


This expression is easily derived from equation (B.5). Note that for this 
calculation the angles 0 and 6 are averages calculated for an observer midway 

1 . 1 1 ^ 

between the two microphones. 

Figure 22 shows the percent error between the exact (At ) and approximate (At ) 

e a 

solutions defined as 
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PERCENT ERROR = |qq 


At0 Atq 


(16) 


The calculations are for the maximum microphone separation used which was 
£ = 0.48 m. The error goes as the square of the microphone separation so that the 
next smaller separation value listed in figure 16 would give an error approximately 
1/4 as large. Even at the large separation, however, the error was very small, 
on the order of a few percent. This justifies the experimental approach described 
earlier . 

Comparison of Measured and Theoretical Refraction Angle Correction - Comparison 
of theory and experiment is given first for the case of an on-axis source location. 
Shown in figure 23 is a comparison between the measured and theoretical angle 
correction as a function of open jet Mach number and acoustic source frequency. 

The source and far field microphone positions correspond to Test Configuration 1. 
Specific coordinates for this geometry are given in figure 23(a). Only two 
curves are shown in each plot to avoid crowding the experimental data points. 

The agreement between theory and experiment is considered good particularly 
at frequencies above 1 kHz where the phase differences are large and measurement 
.accuracy is correspondingly good. At 1 kHz, the small measured phase difference 
introduced some scatter into the experimental data points. To obtain a quantita- 
tive assessment of the refraction correction theory, the absolute difference between 
the theoretical curves and the data can be calculated. Considering the measurement 
angle, 0 = 37.5 , where the correction is significant, the average difference was 

1.25^ at M = 0.1, 0.25^ at M = 0.2, 0.4^ at M = 0.3, and 0.6^ at M = 0.4. The 
largest deviations occurred at low Mach number due to the small measured phase 
difference. 

Figures 24 and 25 show comparisons between the measured and theoretical angle 

correction at two additional downstream locations. The frequency dependence at 

M = 0.4 is shown in part (a) and (b) of each figure while the Mach number dependence 

at a fixed frequency of 5 kHz is shown in part (c) of each figure. This limited 

frequency and Mach number study was conducted to verify the trends obtained in 

figure 23. In all cases, the agreement between theory and experiment is good. 

The absolute difference between theory and data points was, in all cases, less 

than 2° at e = 37.5°. 
m 

Considering now off-axis source locations, figure 26 is a comparison between 
the measured and theoretical angle correction as a function of open-jet Mach number 
and acoustic source frequency. The test conditions are identical to those evalu- 
ated in the on-axis geometry presented in figure 23. This approach permitted 
assessing the capability of the theory to predict changes in source-to-shear layer 
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separation distances. The results for f = 1 kHz are not included due to the large 
scatter in the data. 

The comparisons in figure 26 indicate good agreement at angles close to the 
jet axis. However, a consistent discrepancy occurred at large measurement angles 
for all frequencies and Mach numbers above M=0.1. AtM=0.1, the discrepancy 

was not perceptible since the angle correction is smaller. 

The origin of the discrepancies was traced to the selection of the source- 
to-shear layer separation distance, h, which is an input to both the experimental 
data reduction (see equation 7) and the theoretical analysis. The results in 
figure 26 (as well as the on-axis comparisons in figures 23, 24, and 25) were 
obtained using a vortex sheet located at the nozzle lip-line to represent the 
shear layer. From the mean velocity calibrations shown in figure 10 this vortex 
sheet was also located close to the half velocity point in the shear layer. The 
value of h/R corresponding to a shear layer located at the nozzle lip-line was 
0.55. If, however, the zero thickness shear layer is arbitrarily assumed to be 
located such that h/R = 0.76, agreement between theory and experiment improves 
significantly as shown in figure 27. In figure 27, the theoretical prediction was 
recalculated using hVR = 0.76, in place of h/R = 0,55, Similarly, the experi- 
mental angle corrections were recalculated using h’/R = 0.76 and equation (7). 

As shown in figure 27 the agreement between theory anS experiment is good at all 
angles with the difference not exceeding 2 at 6 =37,5 over the complete 

operating range. Based on the measured mean velocity profiles, the assumption of 
h/R = 0,76 implies the zero thickness shear layer is being placed at approximately 
the quarter velocity point in the shear layer; that is, U = 1/4 at the assumed 

vortex sheet location. 

This result warrants further discussion since a forced agreement between 
theory and experiment does not provide assessment of theory. By comparing figure 26 
for h/R = 0.55 and figure 27 for h/R = 0.76 it can be seen that there is a 
negligible change in the theoretical angle correction. That is, the theory is not 
particularly sensitive to small changes in h/R . The procedure to reduce experi- 
mental data, however, is sensitive to the assumed shear layer location since the 
ray outside the open-jet is traced back to the shear layer crossover point. As the 
assumed shear layer location is changed, the crossover point shifts and the experi- 
mental result is altered. The better agreement in figure 27, therefore, is due to 
a change in the experimental result. This indicates that the present experiment 
and data reduction procedure are poorly designed for assessing theory at small 

h/R . 
o 


To support this conclusion, the sensitivity of experimental results to a 
small change in h (Ah.) was calculated and found to be significantly greater at 
small h/R (0.55) than large h/R (1.44). As discussed subsequently, good agreement 
between theory and experiment was obtained at the larger value of T^/R^* 
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Locating a zero thickness shear layer on an extension of the nozzle lip-line 
for the purpose of reducing experimental data involves two assumptions. First, 
shear layer thickness is not important and second, the nozzle lip-line represents 
the appropriate location to place the vortex sheet. If a zero thickness layer is 
not assumed, however, refraction theory must be applied to trace the ray back 
through the shear layer as part of the data reduction procedure. This then 
introduces the problem of employing theory, to some extent, to reduce data which 
in turn is intended to assess theory. Such a calculation would be of value, however, 
in establishing a consistency check (closure of the problem) . In addition, a 
detailed error analysis of the experimental procedure to determine its sensitivity 
to changes in assumed values of h would further clarify the problem. 

To further investigate the observed disagreement between theory and experiment 
for the off-axis case, additional measurements were performed with the source at 
the same axial location (X/R = 1.33) but displaced in the opposite direction from 
the jet centerline (h/R = 1.44). Here the source and far-field microphone 
positions correspond to Test Configuration 5. Shown in figure 28 is a comparison 
between the measured and theoretical angle correction at selected Mach numbers and 
acoustic source frequencies. The agreement between theory and experiment is 
observed to be good at all operating conditions. 

In summary, data obtained at three on-axis, axial source positions, four 
source frequencies and several Mach nubmers between 0.1 and 0.4 suggest that the 
refraction angle correction is independent of shear layer thickness and frequency 
as predicted previously. Off-axis theory was validated for source position 
having a large value of h/R . Failure to obtain agreement at an off-axis position 
having a small value of h/R is believed associated with assumptions employed in the 
reduction of experimental data. Based on work obtained to date, therefore, the 
refraction angle change correction procedure should be applied as in the present 
study, that is, the zero thickness shear layer should be located on the nozzle 
lip-line for all source positions. 


Amplitude Correction 

Experimental Approach - The primary objective of the present research was to 
experimentally assess the refraction angle correction theory. A lesser effort was 
devoted to experimentally evaluating the amplitude correction theory. Tests were 
conducted using the same two dimensional geometry (fig. 15) employed in the angle 
correction study. A dipole source consisting of a cylinder undergoing vortex 
shedding was placed in the open jet airstream. The resulting far field tone 
directivity was measured at various Mach numbers. Each directivity curve, repres- 
ented by the function SPL (0 , M), was normalized relative to the sound pressure 
level at 0 = 85^. Theoretical curves for the dipole directivity inside the 
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open jet were also generated for each test condition. These curves were identified 

as SPL (0 (0 ) ,M) . Here, the angle 0 inside the airstream was linked to 0 by 
- - c m , 1 c m 

the angle correction theory. 

Finally, the experimental amplitude correction function, AMP, was determined 
from the expression 


AMP(^m,M) = SPL (^c(^m),M)-SPL(0m ,M)+C|(M) 


(17) 


The parameter, was needed to adjust the right side of equation (17) to agree 
with the absolute value of the amplitude correction. 


Ideally would be determined from a direct measurement of the dipole sound 
pressure level at one specific angle outside the airstream and the corresponding 

angle, 0c(0m) > inside the open jet. However, such measurements exceeded the 
exploratory nature of the present study. Instead, the experimental amplitude 
correction function was equated to the theoretical amplitude correction at the 85® 
reference angle. 

The experimentally determined amplitude correction function was then compared 
to the theoretical prediction over a range of measurement angles. It should be 
noted that this approach is only capable of verifying the shape of the amplitude 
correction theory at each Mach number. Absolute sound pressure levels cannot be 
ascertained since was not evaluated by direct measurement. 

Comparison of Measured and Theoretical Amplitude Correction Results - The solid 
curve in figure 29 shows the predicted amplitude correction based on the theoretical 
formulation of the refraction problem presented earlier. The result applies for an 
acoustic source situated at the coordinates defined by Test Configuration 6. 
Satisfactory agreement between theory and experiment exists at measurement angles 
greater than 55 degrees. The poor agreement at angles close to the jet axis was 
typical of the results obtained during the present study. 

The discrepancy was attributed to the dipole sound source employed in the 
experiment rather than the theoretical analysis. It is believed that the three- 
dimensional flow near the ends of the 20 cm long test cylinder (fig. 6) resulted 
in a spanwise varying vortex shedding process. This condition would not simulate 
the two-dimensional vortex shedding process for which the dipole directivity was 
calculated in equation 17. 
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The above experimental approach was unable to verify directly the theoretical 
amplitude correction. The theory can, however, be verified indirectly using the 
theoretical and experimental amplitude correction study recently reported by Ahuja 
et al. (ref. 14). The theory in reference 14 was assessed for M <.0.26 with the 
results confirming the validity of the amplitude correction factor for a source 
situated on the jet axis. Since the analysis developed in the present study is 
almost identical (to be discussed below) to the theory in reference 14, the present 
analysis was indirectly confirmed for a source on the centerline. It still remains 
to verify the present theory at higher Mach numbers and off-axis acoustic source 
positions . 


Comparisons With Other Investigators 

It is worthwhile to clarify the relationship between the present refraction 
theory and that described in the publication of Ahuja (ref. 14). The angle 
correction for an on-axis acoustic source reported by Ahuja is exactly that given 
earlier by Amiet (ref. 6) and Jacques (ref. 7) for a zero thickness shear layer. 

The general analysis in the section titled Theoretical Formulation of the Refraction 
Problem reduces to that given in reference 9 for the special case of an on-axis 
acoustic source. 

The general amplitude correction theory in the present study also assumes a 
zero thickness shear layer. In contrast the on-axis analysis used by Ahuja assumes 
a shear layer thickness which is infinite compared to the acoustic wavelength. The 
difference between these two approaches for an on-axis source was given earlier by 
Amiet (eq. (14) of ref. 9) and is quite small over most of the angular range, as 
shown in figure 30. At the extreme measurement angles (near the zone of silence and 
near the upstream axis) where the difference between the two approaches becomes 
significant, it would be expected that neither result would be accurate since neither 
approach properly accounts for the multiple reflections which occur in these regions. 
Thus, in the range of practical use, the approach of Ahuja gives results which differ 
little from those of Amiet (ref. 6) or the present analysis. Consequently, the 
experimental confirmation of the theory in reference 14 also confirmed indirectly 
the present analysis. The advantage of the present analysis is that it can treat 
the case of an off-axis source. 


Evaluation and Summary 

Refraction of sound waves passing through an open jet circular shear layer 
results in significant wave front angle changes at Mach numbers greater than 0.1. 
Acoustic source directivity patterns are thereby changed and, hence, may alter the 
conclusions drawn from an experiment. The refraction angle change is independent 
of the acoustic source frequency and the axial source location. 


30 


A theory capable of predicting the circular open jet refraction angle changes 
for an arbitrary acoustic source type, source position, and far-field observer 
position was validated experimentally. The lip-line vortex sheet model employed 
in the theory is, therefore, an adequate representation of the finite thickness 
shear layer. Axial variations in the mean shear layer properties and divergence of 
the flow field are not needed in the theoretical formulation. The frequency inde- 
pendence of the analysis was verified. Failure to obtain agreement between theory 
and experiment for an off-axis source at small source-to-shear layer separation 
distance is believed to be associated with assumptions in the data reduction 
procedure. 

The theory for the refraction amplitude changes associated with sound trans- 
mission through the shear layer could not be verified experimentally. The three- 
dimensional unstable vortex shedding process from the circular cylinder noise source 
failed to generate the necessary dipole acoustic source directivity pattern. 
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TURBULENCE SCATTERING EXPERIMENT 


Formulation of the Problem 

The objective of the experiment described in this section was to assess the 
importance of discrete tone frequency scattering and spatial scattering as sound 
propagates through the open jet turbulent shear layer. The first mechanism extracts 
energy from the tone and redistributes it in adjacent frequency bands resulting in 
a broadened spectrum while the second mechanism redirects sound to new angles. Both 
mechanisms result in a change of the discrete tone directivity pattern. 

Measurements obtained during the angle correction experiments indicated that 
frequency scattering occurred in the UTRC open jet. For example, spectrum broadening 
was found for the source and microphone location defined by Test Configuration 1 in 
Table I . Figure 31 shows the result with f = 10 kHz, M = 0.4 and 0^^ = 90^. Here 
the sound was transmitted through the single shear layer between the source and the 
far f ield microphone . 

Unfortunately, the experimental configuration used in the angle correction 
measurements in figure 5 cannot be used to directly assess discrete tone scattering 
phenomena. Both refraction angle and amplitude changes occur across the shear 
layer making it difficult to isolate the tone amplitude changes due to scattering. 
Also, varying the open jet Mach number can change the speaker source strength and 
directivity pattern inside the airstream. Thus, changes in tone amplitude with 
Mach number cannot be attributed, necessarily, to turbulence scattering phenomena. 

To circumvent these problems an alternative test configuration was chosen. 

An acoustic source was located on one side of the open jet (figure 7) transmitting a 
pure tone through both shear layers to an array of microphones on the opposite 
side of the airstream. In this case, the speaker source strength could be assumed 
constant with changes in tunnel Mach number. The tests were limited to assessing 
scattering effects at 90^ to the open jet. This technique had previously been 
employed by Paterson (ref. 20) to investigate turbulence scattering in a small, 
turbulent round jet. 

The first microphone, m^, was situated on the acoustic source centerline and 
measured the zero Mach number sound pressure level. The remaining microphones were 
displaced downstream at a distance £, to account for convection through the open 
jet. Analytically, the displacement relative to the m^ microphone was given by 
the simple Mach number dependent expression £ = 2 RqM. Based on the above test 
configuration, microphone m^ measured the M = 0 acoustic spectrum while m^ measured 
the M = 0.4 spectrum. The experimental approach in figure 7 assumes that microphone 
m^ is situated on the speaker centerline. If this condition is not satisfied, 
the measured sound pressure level will correspond to the intensity on a side 
of the lobe of the speaker directivity pattern. 


32 


The relative importance of any microphone misalignment with the speaker 
centerline is determined by the shape of the speaker directivity pattern. This 
was evaluated by a separate test in which a microphone was traversed on a line 
parallel to the array in figure 7. The speaker centerline was then identified as 
the station at which the intensity was a maximum. The changes at locations off the! 
speaker axis were then tabulated after correcting for the increased distance to ! 
each microphone location. At 5° from the speaker centerline, the following changes 
in the directivity pattern were observed: 

Frequency Change in SPL 

1 dB decrease 

2 dB decrease 

1 dB increase 

2 dB increase 

The sound pressure level increased or decreased monotonically between the speaker! 
centerline and the 5° off axis position. 

The geometrical technique used to locate m^ in figure 7 was capable of aligning 
the microphone within 5^ of the acoustic source centerline. Considering the^ 
small change in the speaker directivity pattern over this range, it can be concluded 
that microphone m^ measured the approximate intensity emitted on the speaker 
centerline. 

Placement of the downstream microphones in figure 7 assumes that the open jet 
can be modeled as a uniform velocity profile bounded by a vortex sheet. In reality, 
this approach fails to account precisely for the streamwise convection of acoutic 
wavefronts as they propagate through the finite thickness open jet shear layer. 
Consequently, the downstream microphones would not be aligned with the acoustic 
ray emitted on the speaker centerline and emerging at 90^ to the open jet. Instead 
these microphones would measure the intensity on the side of the lobe of the 
speaker directivity pattern. 

The relative importance of propagation through the shear layer can be 
estimated by modeling the open jet as a uniform velocity profile bounded by a 
vortex sheet at the 10 percent velocity line in figure 10. At X/Rq = 4 and M = 0.4, 
microphone m 4 in figure 7 would be 7^ off-axis relative to the acoustic ray 
emitted on the speaker centerline. On the other hand, if the open jet diameter 
corresponds to the 90 percent velocity line, then m 4 would be 4^ off-axis. In 
either case the spear directivity pattern is sufficiently broad to ensure that 
microphone m 4 measures approximately the intensity of sound emitted on the 
speaker centerline. Thus, modeling the open jet as a uniform velocity profile 
bounded by the lip line vortex sheet is adequate for locating the downstream 
microphone positions in figure 7. 


5 kHz 
10 kHz 
20 kHz 
25 kHz 
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Admittedly, the dual shear layer transmission in figure 7 creates more severe 
frequency and direction scattering than the single shear layer transmission in 
figure 5, This is verified by the spectrum comparison in figure 31. Consequently 
the results presented in the following section represent trends rather than direct 
calibrations for turbulence scattering through a single shear layer in an open jet. 


Experimental Results 

Figure 32 shows the measured discrete tone spectra as the open jet Mach number 
and source frequency changed. All spectra (with the exception of those in figure 
32(d)) are plotted relative to the tone sound pressure level measured at M = 0. 

This permits monitoring changes in tone amplitude due to both frequency and direction 
scattering. Normalized sound pressure levels were used in figure 32(d) since the 
speaker output amplitude drifted with time. The analyzer bandwidth, BW, increased as 
the acoustic source frequency increased. This was a limitation of the analyzer 
employed in the study. Standard correction for analyzer bandwidth could not be 
made since each spectrum contained both pure tone and narrowband random noise. 

Methods exist for estimating the correction, but such detailed information exceeded 
the exploratory nature of the present experiment. 

Several conclusions can now be arrived at based on figure 32. First, tone 
broadening, at frequencies above 5 kHz, increases as the open jet Mach number 
increases. Second, tone broadening is accentuated as the source frequency increases. 
Finally, the tone amplitude decreases as both frequency and Mach number increase. 

Few open jet experimental studies of acoustic source scattering effects are 
available for comparison. Abuja et al. (ref. 14) have shown that for = 90^ and 
frequencies below 10 kHz, no broadening occurs about the discrete tone peak value. 

This result, which is limited to M _< 0.18, agrees qualitatively with the observations 
in the present study for similar test conditions. Ahuja also found that the discrete 
tone amplitude remained almost constant if f _< 10 kHz and M _< 0.18. Again this 
agrees qualitatively with the results in figure 32. Although a 2.5 dB decrease 
occurred in the present investigation at f = 10 kHz and M = 0.2 this may be due to 
the dual shear layer transmission approach employed here. 

The importance of changes in turbulence length scale was also investigated in 
the present experiment. This was accomplished by locating the source at three axial 
stations given by X/Rq = 1.33, 2.66 and 4.00. Turbulence length scales are known to 
scale with the local shear layer thickness. Based on the shear layer growth rate 
in figure 11 the eddy length scales increased by a factor of 1.7 as the acoustic 
source moved downstream. Figure 33 shows the resulting increase in frequency 
scattering and a decrease in absolute tone sound pressure level at f = 5 kHz, M = 0.4. 
Similar results were observed at higher frequencies. The measurements demonstrate 
that turbulence scattering effects are more pronounced at downstream stations. 


34 


Based on the above trend, shear layer scattering at stations far downstream 
must be assessed accurately if source directivity patterns are required for discrete 
tones. This follows since sound radiating to the far field at angles close to 
the jet axis is originally transmitted through the shear layer at large X/Rq. 
Consider, for example, an acoustic source situated on the open jet axis in Test 
Configuration 1. If M = 0.4 and 0^ = 22.5^ then, based on the angle correction, 
the shear layer crossing point is given as X/Rq = 3.2. Consequently, sound 
radiating at angles close to the jet axis encounters more severe turbulence 
scattering than sound radiated at 90^ to the open jet. 

Ahuja observed that the ratio of shear layer thickness (or eddy size) to 
acoustic wavelength was not large enough for significant scattering to occur. This 
statement must be considered in light of the low Mach number (M <_ 0.18) and the 
axial acoustic source location employed in Ahuja's study. The latter parameter 
was limited to one location given by X/R^ = 0.44. Based on the Goertler sheari 
layer thickness theory (fig. 11), which would be expected to describe the open 
jet shear layer thickness used in reference 14, the eddy length scales at X/Rq = 0.44 
would be small compared to the present experiment. Consequently, scattering effects 
would be weaker. However, scattering effects would be more pronounced with changes 
in axial location in the Ahuja study since the shear layer growth rate is larger. 

At X/Rq = 4.0, both facilities would be expected to demonstrate similar scattering 
effects. 

Although scattering was observed at high frequencies and high Mach numbers in 
the present study, it was not considered an important phenomenon for the angle 
correction experiment. This is because the wave front tracking technique employed 
in the present study is independent of absolute sound pressure level. Admittedly, 
spatial scattering from points adjacent to the shear layer crossing point, X^, in 
figure 15, would contribute to the cross-correlation function. However, the 
contribution would be random since the eddy size and instantaneous velocity is 
random. Long time averages would tend to smooth out these contributions leaving 
only the refraction effects due to the mean flow. Similar arguments can be made 
for the frequency scattering effects. The refraction amplitude correction 
experiment was not susceptible to scattering effects since the dipole source 
frequencies were below 5 kHz. 


Summary 

Discrete tone broadening and absolute sound pressure level changes are signi- 
ficant at frequencies above 10 kHz and open jet Mach numbers greater than 0.2. These 
observations apply if the source is within one diameter of the open jet exit. 
Scattering effects are more pronounced at stations further downstream with scattering 
phenomena appearing at frequencies below 10 kHz. 
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The above observations demonstrate that turbulence scattering can be 
significant in open jet test facilities at Mach numbers approaching M = 0.4, 

If accurate source directivity patterns are required for discrete tone sources, 
the scattering mechanism must be considered. Model propeller and rotor noise 
studies as well as supersonic jet screech noise investigations, therefore, can be 
adversely affected by turbulence scattering phenomena. 
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CONCLUSIONS 


1. Far-field noise directivity patterns measured in open-jet acoustic test facilities 
are significantly altered at test Mach numbers of 0.1 and greater due to sound 
wave refraction by the open- jet shear layer. This is evident from measurements 

of the wavefront angle change associated with sound propagation through the 
shear layer. 

2. For on- axis source locations, the refraction angle change for a circular cross- 
section open jet was well predicted by the zero thickness shear layer theory of 
Amiet over the Mach number range of 0.1 to 0.4 employed in the present experiment. 

3. For the range of shear layer thicknesses considered in the present experiment, 
the refraction angle change was independent of shear layer thickness. This is 
evident from wavefront angle change measurements conducted at axial source posi- 
tions between 0.3 and 2 diameters downstream of the nozzle exit. This independence 
of thickness and shear layer divergence confirms a previous theoretical prediction. 

4. The refraction angle change was independent of frequency over the 1 kHz to 10 
kHz range considered in the present experiment. This independence confirms a 
previous theoretical prediction. 

5. The off-axis source, refraction angle theory developed in the present program 
was confirmed experimentally for the case in which observer and source were 
located in a plane passing through the open-jet axis and the source- to-shear 
layer distance was greater than one open jet radius. 

6. Measurable differences between theory and experiment for the refraction angle 
change occurred with an off-axis source at a source-to-shear layer separation 
distance less than one jet radius. This disagreement, at present, is not 
viewed as a failure of the theoretical angle correction. Rather, it is believed 
due to the sensitivity of the experimental method to small changes in the 
assumed zero thickness shear layer position at small source-to-shear layer 
separation distance. 

7. Theory for the refraction amplitude change associated with sound transmission 
through a shear layer could not be verified due to limitations of the experimental 
technique. Comparison of predicted refraction amplitude changes with previously 
published measurements provides a partial validation of theory for small Mach 
number and on- axis source locations. 

8. Scattering of sound by shear layer turbulence is significant at Mach numbers 
greater than 0.2 and frequencies greater than 5 kHz. Scattering effects increase 
as the source position is moved downstream. 

United Technologies Research Center 

East Hartford, Connecticut 06108 

April 15, 1978 
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APPENDIX A 


Derivation of Shear Layer Correction For An Off-Axis 

Acoustic Source 


Figures 13 and 14 show the basic geometry and the coordinate system used in the 
following derivation. These figures were previously described in the section entitled 
"Theoretical Formulation of the Refraction Problem." 

Refraction Angle Change - The angular change of a ray on passing through the 
shear layer is determined by treating the waves on either side of the shear layer 
as plane waves and matching phase across the shear layer. This planar assumption 
should be valid as long as the wavelength is smaller than the open jet radius. 

The small disturbance pressure field produced by an arbitrary source can be 
written as a Fourier composition of plane waves of the form 


p_ - g-i(0-k2,z) 


(A.l)l 


where 


@s-cut +k5(X+kyy 


and k and represent wavevectors in the x, y, z coordinate system, 
tion 2f equation (A.l) in the wave equation 


(^-Co = 0 


+ IJ — 

Dt " dt ^ ^ dx 


(A. 2) 
Substitu- 

(A.3) 


gives 


k,| = (l-MkJ^- k/-k^ 


(A. 4) 


where the overbar on k indicates normalization by o)/c ; i.e., £ = k C /(o. 

o o 

The normal to a wavefront is found by taking the gradient of the phase of 
equation (A.l) giving 

7 ^ r 'N r r 

A ky 1+ ky j -b k^|k (A. 5) 


n. 


I-Mkw 


38 



Since the wavefronts produced by a source in a stream are spheres with their center at 
the retarded source position, this normal vector points from the retarded source 
position to x^, 2 :^'. The direction of energy propagation (the line 

connecting the present source position to x^, found by adding 

i M to equation (A. 5). Thus, 

A Hg+tM i(M +/3^kx)+t ky+ kkz, 

Hi |A a I “ - V • / 

yi+M^[(i-Mk/-k/] 

The transmitted wave at point x^, written 

r j 0 “* (® (A. 7) 

Substitution in equation (A. 3) but with M = 0 gives 



(A. 8) 


Taking the gradient of the phase of equation (A. 7) gives for both the normal to a 
wavefront above the shear layer and the direction along which energy is 
propagated 


^ I 

Hg = kj( 


- A 
I + 


r ^ 

ky j + 


_ A 
kjjk 


(A. 9) 


Equations (A- 6) and (A *9) give the relation between the ray paths for the 

incident and transmitted rays. Rather than have this relation in terms of the 

parameters and k (which must be equal in equations (A. 6) and (A. 9) in order for 

the phase of the incident and transmitted waves to be equal across the shear 

layer) it is more useful to put the relation in terms of the angles 0 and 4). 

Then, x , y, and z can be written as 
1 1 1 

X, = r, COS 6(. 


y, a sin 9^. cos 
Z| = r, sin 6^ sinXg 


with 


a sin a a (sin - gsin 

cos<^gSin0Q sin 0Q 


(A. 11) 
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I 


and 


/ = + a 

g = I - h/a 


The first equality in equation (A. 11) can be verified from figure 14, and the 
second equality is found from equation (A. 20) , to be derived later. The normal n 
along r^ is then 


ri| = tcos 9^+^ sin 9^ cos + k sin 9q siny^ 


(A. 12) 


Comparisons of equations (A. 6) axKi (A. 12) show that 


or 


cos 6r = 


M + 


^ - k/] 


(A. 13a) 


I 


COS dc 


’'“’/s'' [ y I - M^sin^Sc 


- M 


Also 


ky = 


sin Bq cos Xj. 

^ I - M^sin^0c 


(A. 13b) 


k„ = 


sin 6r sin 


'c sin Xc 


^yi-M^sin^^c 


(A. 13c) 


For the transmitted wave the normal n 2 along r^ can be written 

A ^ A A 

ng = i COS ©I + j sin cos Xj + ^ sin 6 ^ sin 
Comparison with equation (A. 9) gives 

kx = COS 


(A. 14) 


(A. 15a) 
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ky = sin 5jCos 


(A. 15b) 


kz 2 = sine, sin 


By equating the values of E and E given by equations (A. 13) and (A. 15), the 
incident and transmitted angles ar? found to be 


cos e, = , 

V l-M^sin^ft 


COS dp 


sin 0* cos 


t Xf 


sin 0g cos y^ 
s/ |-M^sin^0c 


Equation (A. 16) can also be written 


ton Br = 

/3 COS0, + M 


where 


C ^ = (I - M COS 0,)^ - cos^ 0, 


which agrees with equation 1 of reference 6. Other useful relations are 

C. I 


sin 0c = 


l + M^C,^ 


sin 0c 

I -M^sin^0r 


The angle a can be related to (|) by applying the sine rule to triangle SAB in 
figure 14. Thus, 
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cos y^ = q cos 


(A. 20) 


Using equation (A.,19) , equation (A. 17) can be rewritten 


sin 6^ COS = COS; 


(A. 21) 


The above relations determine completely the angle change needed in the shear 
layer correction. 

Refraction Amplitude Correction - The amplitude correction will now be determined 
by calculating the ray tube divergence. Consider a ray tube beneath the shear layer ^ 
formed by varying 0^ and by amounts d6^ and d(j)^, respectively. This produces a 
ray tube with cross-sectional-area at the shear layer 


dA| = sin d 6^ d<^ 


Above the shear layer the angle variations are 


39. 

’ dd(. 

dA = 11^- d©c + — d<f> 

♦ dSc ^ dA 


(A. 22) 


(A. 23a) 
(A. 23b) 


Note from equation (A. 16) that 96j./3(()^ = 0 since 6^ and are independent variables. 

Therefore, this term does not appear in equation (A. 23a). The partial derivatives 
can be found from equations (A. 16) and (A. 21) to be 


d^c sin0fSin^5^ 


^ ^ gCtSinc^^ 

d<#>c sin sin 9^ d(j>^ 


(A. 25a) 
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(A. 25b) 


da sin 4>r 

_ g LC 

^4>c sin 


d<f>^ cos 7^ cot 0^ sin 0J.- cos 7c cos 0c (A 26) 

sin 7^ sin0|(l-M^sin2g^)V2 


In order to find the area which the ray tube intercepts on the tangent 
plane, note that will remain constant in equation (A. 10) as 6 and (p are 
varied. From equation (A. 10) ^ ^ 

X| _ cot 0c 

Sin 7^ 

4; = cot 7 ^ (A. 27) 


Thus, 



(A.28) 


A variation d0^ gives a contribution to dx , while a subsequent variation 

d(j)^ gives both a dx^^ and dy^ contribution.^ This indicates that the ray tube 

area intercepted by the tangent plane is a parallelogram; its area dA^ is the 

do contribution to dx times the d(J> contribution to dy , or 
c 1 c 1 


dA - ^ 

^ sin 7 ^ sin 7 ^ sin 0c 


(A. 29) 


The cross-sectional area dA^ just above the shear layer is then, by a similar 
argument 


dAj = dA 2 sin 7 ^ sin 0^ 


(A. 30) 
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Equation (A. 30) can also be derived by setting r2 = 0 in equations (A. 34) to be 
derived later. 


The rate of ray tube divergence changes on passing through the shear 
layer, and the ray tube cross-sectional area must now be calculated for a 
point at an arbitrary distance above the shear layer. The equations for a 
point x^, y^, are 

Xg = X| + fg cos 8 ^ = r cos ©rn 
Yg = y, + fg Sin 8 ^ cos 


with = Xg + Yg + z| 


Zg = Z| + Tg sin 8^s\n 


By varying first 0 and then (j) , the point x , y^, Zg will trace out two 
sides of an area. By taking the cross product of^the vectors defining these 
two sides and then taking the dot product with the unit vector n^ given by 
equation (A. 14), the ray tube cross-sectional area will be determined. In 
taking the derivatives of equation (A. 31) , V 2 will be taken as constant, as the 
variation of r„ with 6 and 6 is irrelevant in the final area result. Thus, 

2 c c 


d8c 

dzg 


dx 


2 - ~''l cin ^ 

d8r sin 6c ^ ^ d8p. 
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I . « 

Tg ^cos 8 ^ cos y^ — Sin sin 

I « 

Tg ^cos siny^ + sin 8 ^ cos 



(A. 32) 


where is constant and the variation of and with 0^ was found using 
equation (A. 28). Also, 

^=-r, cos 0c cot 

sin 0c d(f>. 

~n~ ^2 sin 0t sin y, - — . , 

' <? n ^ ^ ^ (A 


di d<p 

— = r, sin 9, cos y, 


The cross-sectional area dA^ is then 


/a ^^2 A 

^ dO^ ^ 


A ^V2 , A ^^2 




l^=('aT 


a -I. i £- + b =1- 

a* I 34,^ a+. 


dA, s Og -|dc, X dCg) 


Thus, 


rg sin \ /. ^ ’’2 sin yc sin^ d<p^ 


J £2 sin0c ^ 

i' sin0t dej V r, 


siny^ sin0c 


(A. 35) 
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Equation (A. 35) gives the ray tube divergence behavior above the shear layer. 
This equation allows one to calculate the sound level just outside the shear 
layer in terms of the far-field sound. 

The incident pressure just inside the shear layer is now needed, and will 
be calculated from the shear layer transmission coefficient. In crossing the 
shear layer, the wavefront can be treated locally as being plane. Then the 
velocity potentials for the incident, reflected and transmitted waves can be 
written 


* = e~' ■•■''21^) 

i 

= R e*' (A. 36) 

(j> = T. e"'^® ''22^) 

T (p 

where 0 is given by eqqation (A. 2) , k 2 by equation (A. 4) and k_ by equation (A. 8). 

1 Z 2 

The pressure and fluid displacement will be matched across the shear 
layer. The pressure is related to the velocity potential by 


P = 



(A. 37) 


Matching the pressure across the shear layer at Zj^ gives! 



iZ|kz, 


+ Re 


z,k 


z,].. 




(A. 38) 


To match the fluid displacement across the shear layer, the interface can be 
considered to be rippled by the acoustic waves, the ripple moving in the x 
direction with velocity k /U . If the observer moves in the x direction with 
the ripple, the mean flow veiocity outside the shear layer will be while 

that inside will be U -o)/k . Denoting the perturbation z velocity by w and 
equating flow slopes gives 


♦ I - Mkx 


(A. 39) 


46 


Then, using the fact that w = 9 $/az and using equations (A. 36), equation (A. 38) 
becomes 




(A. 40) 


Combining equations (A. 38) and (A. 40), 


[l + (l-Mkxf k^gAzi] = 2 (l-Mkx) e'^i^''22'k2i), 


(A. 41) 


T^is the transmission coefficient for the potential. The transmission coefficient j 
for pressure is found by noting from equation (A. 37) that 


- 

Pf 



. Tol 1 

Pi 


(l-Mkx)^i 

l-Mkx ! 


Thus, 


Tn = 


I + (l - Mkx)^ k^g/kz, 


From equations (A. 13c), (A. 15c), and (A. 19) 


k^2 sin sin 

kzi ■ CfSinxc 


(A. 42) 


(A. 43) 


(A. 44) 


The final result for the amplitude correction can now be given. Equation 
(A. 35) gives the ratio of the square of the pressure just outside the shear layer 
to that at the observer. Equations (A. 43) and (A. 44) give the ratio of pressure 
just outside to the incident pressure just inside the shear layer. Thus, the 
incident pressure just inside the shear layer can be calculated from that 
measured at the observer position. By extrapolating this calculated pressure 
by the factor r^/r, the corrected pressure at a distance r from the source can 
be calculated. The final result for the correction factor is then 
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Pc 

j 

(iL + ’’2 \ 

/r, 1 rg sinyc sin0^ \ 


1 

v '■ sin20jSin0c/ 

\ r r sinyj sin 0c d4>^ I 


(A,45) 


r3/2 

/siny^d^^/^ 

2 sin 0c v sin 0^ sin 


, r. x2 sin 0.sin y. 

+ (l-MCOS0f) i ^ 

CfSinyp 


r, t 2 ( sin 




sin 0c 
r 3/2 


t 

sin 0jsin yj 
sin yca<^^/d0c J 


where 


9<l>t 


9<|)t 


are given by equations (A. 25a) and (A. 26), respectively. 
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APPENDIX B 


Wavefront Angle Change Calculation 
According to equation (5) derived earlier 


a't“ At = a'-a = 





( 5 ) 


Using the trigonometric identity for the sum of two angles and expanding in a 
Fourier series gives 


a't-At= 


Mo' 

cos^, (l — 2“ ) sinyu, 


-( 

Vo 3 


)] 


(B.l) 


Since was shown to be small (y - .016 radians), higher order terms may be 
neglected leaving 


a't - At 



I . Mo(cosMi-i) 
' sinyLt, 


(B.2) 


or 


sin/Lt, 



(B.3) 


The maximum contribution from the second term in equation (B . 3) occurs 
when y^ = tt/ 2. Then 


A 


I 

T 



(B.4) 


Here y provides a negligible contribution to the phase difference. Hence, 
an adequate approximation for the phase difference is 

a't-At = ’^sinyu., (B.5) 

This expression is valid to within .016 cycles. Finally, 



(B.6) 
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APPENDIX e 


Shear Layer Angle Correction Calculation 

Figure 35 depicts the ray path as the acoustic wave is transmitted through 
the shear layer. The point, X^,at which the sound emerges from the airstream 
is given by 



(C.l) 


But 


X 2 = r sin (90 


(C.2) 


and 


X 3 =(Y-h) ton 


(C.3) 


Thus 


Xq = r sin (90-fiJ +{Y-h) tan(/i, -90 + 


or 


Xq = r cos 9ni -(r sin 8^ -h) cot (^, + dm) 


(C.5) 


Finally the original propagation angle inside the flow is given by the expression 




(C.6) 


APPENDIX D 


Microphone - Acoustic Driver Cross-Correlation Technique 


Facility background noise contributed significantly to the two microphone 
cross-correlation function at Mach numbers greater than 0.1. This effect can be 
demonstrated analytically. The net pressure, P, sensed by microphone m 2 in figure 
35 is 


Pglt) = Bg COS [wt - 27t(A + Ag)] + P2(t) 


(D.l) 


Here B 2 represents the frequency dependent discrete tone amplitude while P 2 
denotes the facility random background noise. The time dependence is expressed 
relative to the signal generator which drives the acoustic source in figure 36 1 
The parameter A represents the total phase lag between the signal generator 
output and the microphone input. Here 


A = Aq+ As + Ap 


(D.2) 


where Ap is the propagation phase delay between the speaker and the microphone. 
A similar expression exists for the signal at m^^ except that the pure tone is 
out of phase by A’ relative to m 2 : 


P, (t) = B, cos [(fcirt - 27 t ( A + Al + A,)] + P,'(t) (D . 3 ) 

The resulting cross-correlation function S, is then 

s, 2 (t) = B, B'g c^cos [wt- 27 r(A +A' +A,)] • 

COS j^*nt-27r{A +Ag)j^ 


The second term in equation (D.5) corresponds to the random noise contribution which 
dominates the correlation trace in figure 18. 


51 


An alternate approach for extracting only the periodic component is to 
separately cross-correlate each microphone signal with the acoustic driver signal. 
This technique, described by Schlinker in reference 21, was applied to the 
measurement of standing waves in a duct using flush mounted wall microphones. The 
microphones sensed both the acoustic pressure fluctuations and the random turbulent 
boundary layer pressure fluctuations. 

The method is illustrated schematically in figure 35.1 The output voltage, 

V 2 , from m 2 is cross-correlated with the signal generator voltage, Vg(t) = Vq coswt, 
resulting in 


$2 gW = Gg V q Bg <(cOS [wt - 27T (A + Ag)] cosftol + t]^ 

r (D • 6) 

Note that the microphone system sensitivity, G 2 , has been used to link the acoustic 
pressure and the voltage output. Replacing the trigonometric terms by the complex 
notation for the cosine function and noting that the time average of the periodic 
terms is zero gives 


= "2 ®2 BgVoCOS [wt “ 2 ir (A + Aj)] 


(D.7) 


The measured cross-correlation shown in figure 20 can be represented by a peak 
voltage amplitude H 2 and a phase between Vg and P 2 


* H 2 COS jwT - 27T (n2+'*72^j 


(D.8) 


where n is an integer. 

Comparing like terms in equations (D.7) and (D.8) provides 

’?2=A+A2-n2 0<^2^' (D.9) 

The value for TI 2 is obtained from the microphone - signal generator cross-correlation 
curve in figure 20. The time delay, T 2 , to the first maximum is used to solve 
for ri 2 * Notice that the location of the first maximum was defined as the midpoint 
between two cross-over points in the correlation function. By using two cross-over 
points, irregularities in the correlation function trace were averaged. 
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Figure 20 represents the same operating condition as in figure 18 .1 The 
background noise which dominated in the first figure has been eliminated from 
the second figure. To verify that the background noise did not influence the 
results, the cross-correlation function for microphone m 2 was measured at each 
test condition with the acoustic source turned off. 

Cross-correlating and Vg^.^it can be shown that 

A' = (77,+n2)-A-A, ,0<7j, < co 

Here the number of integer cycles phase delay is referenced to m 2 * Hence, r \2 
can exceed one cycle. 

Finally, 


A't-At = [A+(A,-A2)]-At 

Substituting equation (D.IO) and (D.9) for A’ and n 2 gives 


A't-Aj =(7^1-772) -At 


(D.IO) 

- Til 

(D.ll) 

(D.12) 
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TABLE I - REFRACTION ANGLE CORRECTION 
TEST CONDITIONS 


Test 

Geometry 

Number 

Axial 

Source 

X 

Ro 

Radial 

Source 

h 

Ro 

Microphone 

Position 

h 

r 

Open Jet Mach 

Number 

M 

Source 
Frequency 
f , kHz 

1 

1.33 

1 


0.1, 0.2, 0.3, 0.4 

1, 2.5, 5, 10 

2 

2.66 

1 


.2 

5 





.4 

1, 2.5, 5, 1C 

3 

4.0 

1 

.25 

.2 

5 





.4 

1, 2.5, 5, 10 

4 

1.33 

.55 

,14 

.1, .2, .3, .4 

1, 2.5, 5, 10 

5 

1.33 

i 

1.44 

i 

.36 

.4 

5, 10 


TABLE II - REFRACTION AMPLITUDE CORRECTION 
TEST CONDITIONS 


Test 

Axial 

Radial 

Microphone 

Open Jet Mach 

Vortex 

Geometry 

Source 

Source 

Position 

Number 

Shedding 

Number 

X 

h 

h 

M 

f, kHz 


^o 

^ „ 

r 



6 

1.33 

1 

0.25 

0.22, 0.42 

1.1, 2.1 

7 

1.33 

.55 

.14 

.22, .42 

i 

1.1, 2.1 


56 



















78 04 32 


58 


Figure 2 — Anechoic Chamber Test Arrangement for Refraction Angle Correction Experiment 










INTENSITY, dB 



Figure 3— Intensity— Distance Variation for Acoustic Source with Discrete Frequency Excitation 
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Figure 4— Acoustic Source and Far Field Microphone Arrangement for Refraction 
Angle Correction Experiment 
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78 - 04 - 32-1 







Figure 5— Plan View of Open Jet, Acoustic Source and Far -Field Microphone 
Arrangement for Refraction Angle Correction Experiment' 


78 - 04 - 32-4 
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VIEW LOOKING UPSTREAM 


Figure 6-Schematic of Dipole Source and Microphone Arrangement for Refraction 
Amplitude Correction Experiment ! 
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Figure 7— Plan View of Open Jet, Acoustic Source and Far-Field Microphone: 
Arrangement for Discrete Tone Scattering Experiment- 




NORMALIZED VELOCITY, U/U, 



Figure 8 — Variation of Normalized Axial Mean Velocity Profiles with Axial Location 
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0.4 0.6 0.8 1.0 1.2 1.4 1.6 

NORMALIZED RADIUS, R/Rq 


(b) X/Rq= 1.33 

Figure 8— Continued 
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Figure 8 — Continued 
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(e)X/Ro = 3.33 


Figure;^— Continued 
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1.6i 



NORMALIZED AXIAL COORDINATE, X/Rq 


Figure 10 — Axial Variation of Average Shear Layer Characteristics 


78 - 04 - 65-2 , 

I 


71 






NORMALIZED VELOCITY, U/Uq 


1.2 



SIMILARITY PARAMETER, 7? 


Figure!12 — Similarity of Velocity Profiles at Various Axial Locations for Mach Number, M = 0.4^ 
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Figurel4— Cross— Sectional Geometry for Theoretical Model 
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Figure 16 - Model of Wavefront Angle Changes at Far-Field Microphone Location 
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CROSS-COR R E LATION 
RELATIVE AMPLITUDE 



r= 0 


Figure'17 — Two Microphone Cross— Correlation Function at m( 


RELATIVE AMPLITUDE 



Figure 18 — Two Microphone Cross— Correlation Rinction at Mach Number, M = 0.4, 
Measurement Angle, 0^1 = 37.5°, Test Configuration 1.. 
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RELATIVE SPL 



3 ' 4 51 6 ! 7 

FREQUENCY, f, kHz 


Figure 19 — Typical Far-^Field Microphone Spectrum at Mach Number, M = 0.4, 
Measurement Angle, 6 ^ = 37.5^, Test Configuration 1.| 
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RELATIVE AMPLITUDE 
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CORRECTED ANGLE,0c,DEG 


120 



Figure 23^ Comparison of Measured and Theoretical Refraction Angle Correction as a Function of Mach ; 

Number and Frequency. Test Configuration l.j 
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MEASUREMENT ANGLE, 0^^, DEG 


(b) f = 10 kHz i 

Figure 23— Continued 
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MEASUREMENT ANGLE, d^, DEG 
(c) f = 5 kHz 

Figure 23— Continued 
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MEASUREMENT ANGLE, 0 ^, DEG 
(d} f = 5kHz 

Figure 23 — Continued 
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Figure 23— Continued 
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120 



MEASUREMENT ANGLE, 0^, DEG 
(h)f=1kHz 

Figure 23— Concluded 
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CORRECTED ANGLE, 0c, DEG 







MEASUREMENT ANGLE, 6 ^, DEG 
(b)f= 2.5, 1 kHz! 

Figure 24 — Continued 
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CORRECTED ANGLE, dc, DEG 
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(a) f= 10, 5 kHz I 

Figure- 25 —Comparison of Experimental and Theoretical Refraction Angle Correction 
as a jFunction of Mach Number and Frequency. Test Configuration 3^ \ 
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MEASUREMENT ANGLE, DEG 
(b)f= 2.5, 1 kHz 


Figure 25> — Continued 
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CORRECTED ANGLE, 0c. DEG 



0 20 40 60 80 100 120 


MEASUREMENT ANGLE, 0c, DEG 
(c) f = 5kHz 

Figure 25 — Concluded 
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TEST CONFIGURATION 4: 


X/RQ-1.33i 


h/RQ=0.55 


h/r=0.14 

THEORY, M = 0.2 


MEASUREMENT ANGLE, DEG, 

(a) f= 10kHz' 

Figure! 26 - Comparison of Measured and Theoretical Refraction Angle Correction 
as a Function of Mach Number and Source Frequency. Test Configuration 4. 
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MEASUREMENT ANGLE, 0ni. DEG 


(b)f= 10 kHz 

Figure 26— Continued 
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MEASURED ANGLE, DEG 


{d)f = 5kHz 

Figure 26— Continued 
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MEASUREMENT ANGLE,0m, DEG 
(e) f=2.5|kHz 

Figure26 -Continued 
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CORRECTED ANGLE, d^, DEG 
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CORRECTED ANGLE, dc, DEG 



MEASUREMENT ANGLE, 0^, DEG 
(a)f = lO.kHz 

Figure 27 -Comparison of Measured and Theoretical Refraction Angle Correction as a Functionj 
of Mach Number and Source Frequency Using Modified Source to Shear Layer Separation: 

Distance. Test Configuration 4. 
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MEASUREMENT ANGLE, 6 ^, DEG 


(b)f = lOkHz 

Figure 27— Continued 
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MEASUREMENT ANGLE, 0^, DEG 
(c) f=5 kHz 

figure 27i-Continued, 
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MEASUREMENT ANGLE, 0^, DEG 


(d)f=5 kHz! 


Figure 27l- Continued 
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MEASUREMENT ANGLE, 0rn, DEG 


{f)f= 2.5 khz 1 

Figure 27— Concluded 
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CORRECTED ANGLE, 0c, DEG 



Pigure!28 — Comparison of Experimental and Theoretical Refraction Angle Correction 
as a Function of Mach Number and Frequency. Test Configuration 5. 
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TEST CONFIGURATION 6: 






• DIFFERENCE BETWEEN THEORY PRESENTED BY 
AMIET (REFERENCE9) AND AHUJA (REFERENCE14) 



[ifcm 




Figure 31 —Effect of Single and Dual Shear Layer Transmission on Discrete Tone Frequency 
Spectrum at Mach Number, M = 0.4, Frequency, f = 10 kHz, andlAxial Source 
Location X/Rq=1.33 
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RELATIVE AMPLITUDE, dB 



F igu re’ 33 h Variation of Discrete Tone Spectrum withlAxial Source Location for Mach Number,| 
M = 0.4, andlFrequency, f = 5 k Hz | 
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v'z'r V 2l<t- ({5!+ A’ 2 )/^) 
V,= V^(t-('6 + A+Ai)/f) 

A“'Aa+ A5 + A^p 


Figure 35^Schematic Showing Phase Lag Between Signal Generator 'andj Far— Field Microphone Output | 
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